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We investigate localization of entanglement of multipartite mixed Gaussian states into a pair of 
modes by local Gaussian measurements on the remaining modes and classical communication. We 
provide a detailed proof that for arbitrary symmetric Gaussian state maximum entanglement can 
be localized by homodyne detection of either amplitude or phase quadrature on each mode. We 
then consider arbitrary mixed three-mode Gaussian states and show that the optimal Gaussian 
measurement on one mode yielding maximum entanglement among the other two modes can be 
determined by calculating roots of a high-order polynomial. Finally, we discuss localization of 
entanglement with single-photon detection. 
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I. INTRODUCTION 

Continuous quantum variables have emerged during re- 
cent years as a very promising platform for quantum com- 
munication and information processing [lj, |2[ . Using co- 
herent and squeezed states of light instead of single pho- 
tons for encoding and processing quantum information 
possesses distinct advantages. The entangled squeezed 
states can be generated deterministically in optical para- 
metric amplifiers and can be detected with very high effi- 
ciency using balanced homodyne detectors. This should 
be contrasted with probabilistic generation of entangled 
photon pairs in the process of spontaneous parametric 
downconversion and comparatively low efficiency of the 
single-photon detectors. The favourable scaling of con- 
tinuous variable systems greatly facilitates preparation of 
complex multimode entangled Gaussian squeezed states. 
Various three- and four-mode entangled squeezed states 
were demonstrated experimentally [H, 0, @] and ex- 
ploited in elementary quantum teleportation networks 
0, HI ■ Properties of specific multimode entangled Gaus- 
sian states such as symmetric states [lol. Till or 
cluster states [IH, [3 [H, EH were theoretically analyzed 
in detail. 

The multimode entangled states could form a backbone 
of an advanced quantum communication network where 
each node contains a part of the entangled system. An 
interesting and important question is how much bipartite 
entanglement between two nodes of the network can be 
localized on average by performing local measurements in 
each node and exchanging the measurement results via 
classical communication. This so-called localizable en- 
tanglement [13] has been originally introduced and stud- 
ied for discrete variable systems, in particular quantum 
spin chains. Recently, we have generalized the notion of 
localizable entanglement to continuous-variable systems 
[Tc| and investigated the localization of entanglement of 
multimode Gaussian states via local Gaussian measure- 
ments. We have proved that for pure states it is always 
optimal to perform local homodyne detections, i.e. pro- 
jections onto infinitely squeezed states. Moreover, we 
have also shown that non-Gaussian measurements such 



as photon counting allow us to localize more entangle- 
ment in certain cases. 

In this paper we investigate localization of entangle- 
ment for mixed multimode Gaussian states. We con- 
sider entanglement localization by Gaussian measure- 
ments that preserve the Gaussian form of the resulting 
two-mode bi-partite state. We quantify the entangle- 
ment by logarithmic negativity which is a computable 
entanglement measure for Gaussian states We first 
present a detailed proof that homodyne detection is op- 
timal for localization of entanglement of mixed symmet- 
ric Gaussian states [l(| HH, 0, H3|- We then consider 
arbitrary three-mode mixed states and, somewhat sur- 
prisingly, find that in this case maximum entanglement 
among two modes can be sometimes localized by eight- 
port homodyne detection, i.e. projection onto coherent 
states, on the third mode. We thus demonstrate that the 
homodyne detection is not always an optimal strategy 
if we deal with general mixed Gaussian states. Finally, 
we go beyond the realm of Gaussian measurements and 
analyze localization of entanglement with realistic sin- 
gle photon detectors that can only distinguish between 
presence or absence of photons in the light beam but 
cannot resolve the number of photons in the beam. We 
show that in principle even such low-efficiency realistic 
single-photon detectors may be sometimes advantageous 
for entanglement localization compared to Gaussian mea- 
surements. 

The paper is structured as follows. In Sec. II we pro- 
vide a brief summary of the main properties of Gaussian 
states. Then in Sec. Ill we define the Gaussian localiz- 
able entanglement and discuss some of its features. Sec- 
tion IV contains a detailed analytical proof of optimality 
of local homodyne detection for localization of entangle- 
ment of multimode generally mixed symmetric Gaussian 
states. In Sec. V we study the localization of entangle- 
ment for generic three-mode mixed Gaussian states. Lo- 
calization of entanglement by single-photon detection is 
addressed in Sec. VI. Finally, the conclusions are drawn 
in Sec. VII. 
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II. GAUSSIAN STATES AND MEASUREMENTS 

We deal with quantum systems with infinite- 
dimensional Hilbert space such as modes of the electro- 
magnetic field. An TV-mode system can be described 
by the quadrature operators xj,pj, j = 1,...,N sat- 
isfying the canonical commutation relations [xj, Pk] = 
iSjk- Arranging the operators into the vector R — 
(xt,Pi, . . . ,xn,Pn) T the commutation relations can be 
compactly expressed as [Rj,Rk] = i^ljk, where fl = 



-f =1 J and 



J 



1 

-1 



(1) 



States of the TV-mode system can be represented in the 
2.ZV-dimensional real phase space by the Wigner func- 
tion [2l|. Among all possible states, the states with 
Gaussian-shaped Wigner function are particularly im- 
portant. These so-called Gaussian states can be easily 
prepared experimentally using squeezers and passive lin- 
ear optics and can be efficiently handled theoretically. 
Each Gaussian state p is fully characterized by the first 
moments of the quadrature operators (i?;) = Tr(piij), 
i = 1, . . . , 2N and by the covariance matrix (CM) j with 
elements 7^ = ({Ai? i; ARj}), where ARi — R L — (R4) 
and {A, B} = AB + BA. For any CM 7 there is a sym- 
plectic matrix S, i.e. a 27V x 2N real matrix satisfying the 
condition S£IS T — il, that brings the CM to the normal 
form :23] SjS T = diag(z/i, Vi, . . . , vn, Vn)- The quan- 
tities Vi > 1, i — 1, . . . , N are the so called symplectic 
eigenvalues and can be calculated from the eigenvalues of 
the matrix Qj that read as {±iv\, . . . , ±wjv} [191 ]. 

Gaussian states of two or more modes can be entan- 
gled. In the case of just two modes A and B the state 
Pab is entangled if and only if its partial transposition 
with respect to mode B p T ^ B is not a positive-semidefinite 
matrix [H, [H, El, Ef|. On the level of CMs the par- 
tial transposition operation is represented by the diag- 
onal matrix A = diag(l, 1, 1, —1) [25[ and the CM jab 
of the state pab is transformed by matrix congruence 
jab = h-lAB^- 7 '• Let /ii,/i2 denote the symplectic eigen- 
values of tab and we assume without loss of any general- 
ity that P2 < The Gaussian state pab is entangled if 
and only if jab violates the generalized Heisenberg un- 
certainty relation 7 + id > [25| or equivalently if and 
only if fi2 < 1- It is convenient to express the CM jab 
in a block form with respect to A\B splitting, 



A C 

1AB =\ C T B 



(2) 



The symplectic eigenvalues of the matrix jab can be 
calculated using the formula [19] 



Mi, 2 



'<5±V£ 2 -4A 



(3) 



where S = detA + detB — 2detC and A = detjAB are 
the so-called symplectic invariants. The entanglement of 



a bipartite state pab can be quantified by the logarith- 
mic negativity Ej^\pab] — log 2 The logarithmic 
negativity has the important advantage that it can be 
computed analytically for an arbitrary mixed two-mode 
Gaussian state and we have 



Ejst[pab] = max(0, - log 2 p 2 ), 



(4) 



which is a monotonously decreasing function of the lower 
symplectic eigenvalue pi of the matrix jab |20( . 

Arbitrary quantum measurement can be described by 
positive operator valued measure (POVM) whose ele- 
ments IP, are positive semidefinite operators satisfying 
the completeness condition . IIj = 1 , where 1 de- 
notes the identity operator. Gaussian measurement on 
some mode C is any measurement whose POVM elements 
IF, are represented by Gaussian Wigner functions. This 
means that any operator IIj is proportional to a den- 
sity matrix of some (generally mixed) Gaussian state. In 
practice, any Gaussian measurement on optical modes 
can be performed using auxiliary modes prepared in vac- 
uum states, passive and active linear optics (beam split- 
ters, phase shifters and squeezers) and balanced homo- 
dyne detection on each mode. The corresponding contin- 
uous Gaussian POVM characterizing the measurement 
on C reads 



U c (a) = -D c {a)U Q c D ] c (a). 



(5) 



Here Dc{a) = exp(aa c 
ment operator and ac z 



a*ac) denotes the displace- 
-j^{xc + Wc) is the annihi- 
lation operator. ILj is a density matrix of a Gaussian 
state with covariance matrix Jq and zero mean values 
of the quadratures, and the displacement a is determined 
by the outcomes of the homodyne detectors. Note that 
Jq is fixed by the structure of the measurement setup 
and does not depend on the measurement outcomes. The 
displacement operation is an irreducible representation of 
the Weyl-Heisenberg group and Schur's lemma together 
with the normalization Tr[ILj] = 1 implies that 



H c (a)d 2 a = 1. (6) 



This ensures the completness of the POVM ([5]). 



III. GAUSSIAN LOCALIZABLE 
ENTANGLEMENT 

We consider TV-mode Gaussian state Pabc shared 
among N parties A, B, Cj, j = 1, . . . , N — 2, with each 
party possessing a single mode. Parties Cj attempt to 
increase the entanglement between A and B by perform- 
ing local Gaussian measurements and communicating the 
measurement outcomes to A and B. We are interested 
in the maximum amount of entanglement that can be lo- 
calized on average in the modes A and B in this way. We 
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will use the easily computable logarithmic negativity to 
quantify the entanglement of mixed Gaussian states. 

The elements of the total POVM describing measure- 
ment on all modes Cj can be written as product of the 
single-site elements, 

N-2 

no(«) = ®n Cj (a J 0, 

where a. = (<xi, • • • , arjv-2)> The density matrix of the 
Gaussian state of modes A and B conditional on the par- 
ticular measurement outcome a can be expressed as 



o-ab((x) = 



Tr c [t AB ®n c (a) pabc] 
Tr [l AB ® He (a) Pabc] 



where Trc denotes partial trace over the modes Cj . The 
entanglement of the Gaussian state 17,43 depends only 
on its covariance matrix and not on its coherent dis- 
placement, because mean values of all quadratures can 
be set to zero by means of local unitary displacement 
operations that do not change the entanglement. The 
covariance matrix of <tab(cx) depends only on the covari- 
ance matrix of the initial A^-mode state "/abc and on the 
covariance matrices 7^ characterizing the Gaussian mea- 
surements on each mode Cj. This means that entangle- 
ment of ctab(cx-) does not depend on a which determines 
only the coherent displacement of the conditionally pre- 
pared state. We have that Ej\f \<jab = E^[(jab{0)}, 
Vex. Moreover, any projection onto mixed Gaussian state 
can be expressed as a Gaussian mixture of projections 
onto pure Gaussian states. Clearly, such incoherent mix- 
ing can only reduce the localizable entanglement. It 
follows that in order to determine the maximum Gaus- 
sian localizable entanglement among modes A and B, 
El.gIpabc], it suffices to optimize the logarithmic neg- 
ativity Bjv" [cab (0)] over all local projections onto pure 
single-mode squeezed vacuum states 



3.1 ect: 

m 



E l .g[pabc] = max{£^[o-As(0)]} , (7) 

where 7^ + U > and det(7^) = 1. The definition of 
logarithmic negativity Q implies that the maximization 
([7]) amounts to minimization of the smaller symplectic 
eigenvalue pi of the covariance matrix of the partially 
transposed state <J^ B (0) over all local projections of Cj 
onto pure Gaussian squeezed vacuum states. 

As a side remark, we note that the local projec- 
tions onto pure Gaussian states are optimal among a 
wider class of POVMs whose elements can be writ- 
ten as a product of single-site operators, lie (z) — 
^fj^TlCjizj)- Here 2 is a multi- index labeling the 
POVM elements and the complctncss of the POVM dic- 
tates that J z Hc(z)dz = 1. Each element TlcAzj) is 
Gaussian but the overall structure of the POVM can be 
otherwise arbitrary. We denote by ctab{z) the normal- 
ized conditionally prepared state of A and B. It imme- 



diately follows that 

E l ,g< max E^[a AB (z)} = E^[a AB (z )}. (8) 

where the maximization has to be carried out over all 
Gaussian POVM elements Uj(zj) and z° labels the 
POVM element which maximizes the entanglement be- 
tween A and B. Using Hj(zj) as a seed element for the 
construction of local Gaussian POVM {5} on mode Cj we 
construct a local Gaussian measurement which achieves 
the localizable entanglement En[<jab(z )] and thereby 
maximizes El.q. 

Recently, the problem of localizable entanglement was 
studied for A-mode pure Gaussian states for which it was 
shown that maximum entanglement is localized by homo- 
dyne detection [18(. There it was also noted that homo- 
dyne detection is optimal in the case of Af-mode mixed 
fully symmetric Gaussian states. A detailed analytical 
proof of the latter statement as well as its generalization 
to A^-mode bisymmetric states is given in the following 
section. 



IV. SYMMETRIC N-MODE STATES 

The fully symmetric states are invariant under the ex- 
change of any pair of modes [13, [H| • Owing to the sym- 
metry these states have the highly symmetric CM of the 
form: 



Tsym 







(9) 



where (5 and e denote symmetric 2x2 matrices. The 
matrices /3 and e can be diagonalized simultaneously by 
local canonical transformations that do not affect the en- 
tanglement so that without loss of any generalit y w e may 
assume that (3 = diag(6, b) and e = diag(ei, £2) 

We want to find the minimum of the lower symplectic 
eigenvalue of the matrix jab associated with state of 
modes A and B after local Gaussian measurements on 
remaining N — 2 modes Cj . The minimization needs to 
be carried out over all such measurements. This opti- 
mization task can be greatly simplified if we notice that 
the state © can be brought by unitary transformations 
Uab and Uc on modes AB and C1C2 ■ ■ ■ Cjv-2, respec- 
tively, into the block diagonal form with N — 2 blocks 
2 x 2 7b = 7c 2 = ■ • ■ = 7c«-2 = P - e and one 4 x 4 
block 



lACt 



( P + e ^2(N-2)e 
\ y/2(N-2)e f3+(N- 3)e 



(10) 



The transformation Uc can be physically realized by the 
following sequence of A^ — 3 beam splitters 



Uc = -BciC 2 



1 



sin 



VN^2 



B Cl c 3 sin 



1 



4 



x B, 



Cl CjV-2 



v/2 



(11) 



while the transformation Uab can be performed by a sin- 
gle balanced beam splitter, i.e. Uab — Bab{k/4). Here 
Bij{9) denotes a standard beam splitter transformation 
between modes i and j 



a,i — > a, cos # + aj sin ( 



(12) 



where (at), k = i,j denote standard annihilation (cre- 
ation) operators. 

The interference on the array of beam splitters trans- 
forms the iV-mode symmetric state into a product of 
N — 2 uncorrected single-mode states of modes B and 
Cj, j > 1, and a (possibly entangled) two- mode state of 
A and C\. In this way our task boils down to finding 
a Gaussian measurement on a single mode C\ that will 
prepare the mode A in a state that will give after mix- 
ing with a state with covariance matrix 7 s = (3 — e on a 
balanced beam splitter Uab the maximum entanglement. 
Projection of mode C\ onto a Gaussian state with CM 
7cf prepares mode A in the state with CM [U HI] 



i A = [3 + e - 2(N - 2)e[/3 + (N - 3)e + 7 £ 



Mi-l, 



(13) 



Without loss of any generality we can assume 7^ to be 
CM of a pure squeezed state, 



7# = U{9)V{r)U 1 (9), 



(14) 



where 



U(6) 



cosf 



smf 
cos ( 



V(r) 









,-2r 



(15) 

Mixing of the states with CMs 7^ and 75 on a bal- 
anced beam splitter yields a fully symmetric state with 
CM kab with the diagonal and off-diagonal blocks of 
the form u = (7^ + 7 s)/2 and C = (7^ - 7s)/ 2 , re- 
spectively. The symplectic eigenvalues of the ma- 
trix kab corresponding to the partially transposed state 
can be calculated with the help of the formula (O where 
S = 2(detoj — det£) and A = detKAB- Making use of the 
formula for determinant of a sum of two 2x2 matrices 
P and Q, 

det(P + Q) = dct P + dct Q + Tt[PJQ t J t ], (16) 

where J is given in Eq. (fTJ), one finally finds the square 
of the minimum symplectic eigenvalue fi2 to be 



det (P - e) A n 



where 



in 



7c 



{eig [(/J-eJ-^^-c)-']} 



(17) 



(18) 



and eig(A) stands for the set of eigenvalues of the matrix 
A. Now we seek A m i n such that for any admissible 7 ^ it 
holds 



(/3-e) *i A {l3-e) * - A mln l > 0, 



(19) 



and there is 7 ^ for which the matrix on the left-hand 
side of the inequality has some eigenvalue equal to zero. 
Multiplying the inequality (fTl?)) from the left and right by 
the positive definite matrix (/3 — e) 2 , substituting for j' A 
from Eq. (| 1 3|) and multiplying it again from the left and 
right by the matrix e _1 one gets the equivalent inequality 



where 



X = 



x - y- L > 0, 



e- 1 [f3 + e - X min (l3 - e)] , 



(20) 



2(N-2) 
Y = (3 + (N-3)e + yM. 



(21) 



Multiplying now inequality (|20|) from both sides by X 2 
then by (X 5 2)2 and finally again by X _ 2 one arrives 
at the equivalent inequality 



Y -X' 1 >0. 



(22) 



Substituting here from Eq. (|2Tj) . using the decomposi- 
tion (fT4")) for 7^ and taking into account that (3 and e 
are diagonal, one can rewrite the inequality \22\t in the 
following simpler form, 

Z = diag(di, d 2 ) + f7(6»)^(r)J7 T (6») > 0. (23) 

Here diag(di, d 2 ) is a diagonal matrix with entries 



2(iV-2)e? 

^■ = 6+(JV-3) ej - 



6 + - A min (fr- cj)' 



3 = 1,2. 



(24) 



The eigenvalues of the matrix Z are lower bounded as 



min[eig(Z)] > min(di, d 2 ) + e 



-2r 



(25) 



and the bound is saturated by choosing 9 = if d\ > d 2 
and 9 = ir/2 if d\ < d 2 . The expression on the right-hand 
side of Eq. ([23)) is further lower bounded by min(<ii, d 2 ) 
and this ultimate lower bound is achieved in the limit r — > 
00 which reveals that the optimal measurement on Ci 
localizing maximum entanglement between modes A and 
B is the homodyne detection of quadrature p provided 
that d\ > g? 2 and of quadrature x provided that d\ < 
ei 2 . The sought quantity A m i n can be derived from the 
condition min(di, d 2 ) = which gives using Eq. (f2~4"|) after 
some algebra 



Amin 1 



mm 

j=l,2 



2c, 



6+(7V-3) ej 



(26) 



A close inspection of the above chain of equivalent in- 
equalities confirms that A m i n given by formula (|26p is 
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indeed the minimum denned by Eq. ([18)) . Before going 
further let us notice that in the above derivation we as- 
sumed the regularity of the matrix e, i.e. dete 7^ 0. By 
continuity, the resulting formula (f2T)|) holds also when 
dct e = 0. Without loss of any generality we can assume 
that ei > 62- The minimum in Eq. (126|) is then achieved 
for £2 and A m i n = 1 + 2£ 2 /[6 + (N— 3)e2]. This completes 
the calculation of the minimum symplectic eigenvalue [i 2 
given by Eq. (TiTf . The Gaussian localizable entangle- 
ment of the fully symmetric Gaussian state with CM (J9j) 
can be determined by inserting fi 2 into the formula (j4]) 
for the logarithmic negativity. 

Note that in the present case of fully symmetric states 
we in fact did not restrict ourselves to the local mea- 
surements and therefore we found the optimal mea- 
surement among all Gaussian measurements on modes 
C. Assuming ej > £2, the optimal measurement is a 
collective measurement of the global phase quadrature 

1 

PC.global = VN _ 2 [luj=l PCj 

be, however, implemented locally by local homodyne de- 
tections of the pcj quadrature on each mode followed by 
summation of the obtained measurement results. 

A necessary condition for localization of a non-zero 
amount of entanglement is that det e < 0. First recall 
that f3 — e and (3 + (A — l)e are valid covariance matrices, 
which implies that 



This measurement can 



det(/3 - e) > 1, det[/3 + (N - l)e] > 1. 



(27) 



If det e > then either £1 > and £2 > (i) or e\ < and 
£2 < (ii). In the case (i) we have using Eq. ([26]) that 
Amin > 1. Taking into account the first inequality (|27jl 
we immediatelly get that > 1 and El,g = 0. Consider 
now the second case (ii). From the second inequality ([27]) 
we obtain a lower bound £2 > (1 — b 2 )/[b(N — 1)] = £2, m in- 
The formula Q17p can be rewritten as follows, 

(b-e 2 )[b+(A-l)e 2 ] 

b+(N-3)e 2 {b - ei) - (28) 

For a fixed ei, [i\ is an increasing function of £2 in the 
interval [£2, m in, 0]. A lower bound on can be obtained 
by inserting £1 = and £2 = £2,min into Eq. (f28|) which 
yields /x| > (Nb 2 - l)/(26 2 + A-3). Taking into account 
that b > 1 we finally arrive at /j% > 1 hence El,g = 0. 

The present derivation of Gaussian localizable entan- 
glement for symmetric states can be simply extended 
to the so called bisymmetric states [l2[ that are invari- 
ant under the exchange of modes A and B as well as 
under the exchange of any pair out of A — 2 modes 
Ci, C2, . . . , Cn-2- Such states are described by the CM 
of the form 



7bisyn 



( 13 


e 


T 




... r\ 


£ 


P 


T 




... T 




T T 


a 




... £ 










£ : 










'■■ $ 



(29) 



where = diag(6, b),e = diag(£i, £2), a = diag(a, a), £ = 
diag(£i, £2) arid t are 2x2 matrices. By means of a bal- 
anced beam splitter Uab and the array of N — 3 beam 
splitters (TlT|) we can bring the state into the block diag- 
onal form with one 2x2 block 75 = /? — e, A — 3 blocks 
2x2 7c 2 = . . . = 7Civ-2 = a ~ C an d one 4x4 block 



1AC 1 



P + e ^2(A-2)r 
v/2(A-2)r T a + (A-3)£ 



(30) 



Now the proof of optimality goes along exactly the same 
lines as in the case of fully symmetric states with the 
only difference that now d\ and d 2 are eigenvalues of the 
(generally non-diagonal) matrix 

V = a + (N ~ 3)£ — 2(A — 2)r T [/3 + e - A min (/3 - e)]" 1 ^ 

(31) 

It turns immediately out that also in the case of the 
considered bisymmetric A-mode Gaussian states optimal 
measurement strategy on modes C localizing maximum 
entanglement between modes A and B is homodyne de- 
tection. 



V. THREE-MODE MIXED STATES 

We have already seen that for A-mode Gaussian states 
possessing full or specific partial symmetry with respect 
to the exchange of pairs of modes maximum two-mode 
entanglement is localized by homodyne detection. The 
same statement holds also for pure A-mode states [l8j . 
It is therefore natural to ask whether there are also Gaus- 
sian states for which we can localize more entanglement 
by projection onto a finitely squeezed state than if we 
would use homodyne detection. In this section we give 
an affirmative answer to this question by constructing 
a three-mode mixed Gaussian state for which maximum 
entanglement is localized by projection of one of modes 
onto the coherent state. 

We start with a generic three-mode mixed Gaussian 
state with CM ^abc expressed via the 2x2 submatrices 
with respect to the A|J5|C splitting as 



1ABC 



A D E 
D T B F 
E T F T C 



(32) 



Our goal is to find optimal Gaussian measurement 
on mode C localizing maximum entanglement between 
modes A and B. For the sake of comparison, we will 
also determine the performance of the optimal homodyne 
detection. This latter problem amounts to optimizing 
the relative phase between the local oscillator and signal 
beam in the balanced homodyne detector such that the 
measurement localizes maximum entanglement between 
A and B. We start with a generic measurement in sub- 
section V. A whereas the case of homodyne detection will 
be solved in the subsection V.B. 



6 



A. Optimal Gaussian measurement 

If we project mode C onto a generic Gaussian state 
with CM then the CM matrix jab specifying the 
resulting Gaussian state of modes A and B reads [271 I2II ] 



7AB 



)-(^)(c+7^r 1 (^ T ^ T )- 03) 













-J 



Expressing this CM in the 2x2 blocks as 



7AB 



.4 a 



(34) 



the symplectic invariants needed to calculate the sym- 
plectic eigenvalues given in Eq. ([3]) read as 



S = detA + detB - 2det£, A = det7As- 



(35) 



Assuming without loss of generality projection onto a CM 
of a pure state (|14p and the submatrix C of CM (f32|) to 
be in the standard form C — cl we get explicitly using 
Eqs. (fl"6|) . (|33|) and |35|) the invariant <5 in the form 



5 = detA + detB - 2detD + 



(detE - detF) 2 + [c + cosh(2r)]Trx + sinh(2r)[xi sin(20) + X 2 cos(26>)] 

1 + c 2 + 2ccosh(2r) 



(36) 



r 



where xi = X12 + X21, X2 = X22 - X11 and \ = 
2E T J T DJF - E T J T AJE - F T J T BJF. 

The second invariant A = detjAB can be conveniently 
calculated using the formalism of Wigner functions. The 
value of the Wigner function of the initial three-mode 
state at the origin is equal to 



W A bc(0) 



1 



7r 3 Vdet7ABC 



(37) 



After projection of mode C onto the pure Gaussian state 
with CM Jq we get the conditional (unnormalized) two- 
mode Wigner function whose value at the origin reads 



Wab(P) 



1 



2tt 



7rVdet 7 ^c / det(rc + 7 M-i } ' 



(38) 



where Tc is the 2x2 submatrix of the matrix Tabc 
Iabc expressed with respect to the A|B|C splitting, 



BC 



t d v E 
rl r B t f 
r£ r T F r c 



(39) 



mode C characterized by covariance matrix C on a pure 
state with CM 7^, and we obtain 



Pc(0) 



1 



2tt 



TrVdetC 



V / det(C-i+ 7 ^- 1 ) 



(40) 



A formula similar to Eq. (j37|) holds for the normalized 
Wigner function of the conditionally prepared state of 
modes A and B, and we have 



W AB (0) 



1 



=. (41) 

P C (0) 7rVdet7^ B v ; 

Upon combining formulas (|3"8]) , PO)) and (|4T1) one finds 
the invariant A in the form 



det(r c + 7c )det 7 ABc 
det^" 1 +7 c / ~ 1 )detC 



(42) 



The normalization of the Wigner function Wab can be 
determined as a probability of projecting the state of 



Making use of Eq. (Tl7o|) finally leads to the following ex- 
pression for the invariant A: 



J 



det7Asc{l + detr c + cosh(2r)Trr c - sinh(2r) [.91 sin(2fl) + g 2 cos(26»)]} 

l + c 2 + 2ccosh(2r) ' 

I 



(43) 



where gi — Tc,u + Tc,2i and 32 = rc,22 — Ten- For tic invariants |36|) and (I43|) are nontrivial functions of 
a generic three-mode Gaussian state both the symplec- the squeezing r and phase 9. Consequently, minimiza- 
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tion of the lower symplectic eigenvalue 1x2 given by Eq. 
([3]) amounts to the finding of the minimum of the func- 
tion / 



2^1 



VP 



4A of two variables r and 9. 



Extremal equations df jdr = and df /dO — for opti- 
mal r and 9 yield after some algebra the following pair 
of equations: 

A {d 3 5f 6 (djS) (dj A) + (dj A) 2 = 0, j = r, 9, (44) 

where the symbol dj stands for partial derivative with re- 
spect to the variable j. On inserting in Eq. (|44|) for 5 and 
A from Eqs. (f3"6"| and (j4"3"|) and performing the respective 
derivatives one gets the optimal phase 9 and squeezing 
r to be solutions of the set of two coupled polynomial 
equations, 



y 



(l+x 2 )- 



i=0 j=0 
6 6 



_L_gg (ft)ylV = , 



(45) 
(46) 



where x = tan#, y = tanhr and the coefficients {Pe)^ 
and {Pr)ij are complicated functions of elements of the 
CM (l3"2")) that are not specified here. The equation (143)) 
can be solved with respect to the variable y as 



2/1,2 



-si ± 



4s s 2 



2s 2 



(47) 



where Sj, j — 0,1,2 are sixth-order polynomials in the 
variable x. Substituting now from Eq. (I47p for y into 
equation (I46[) we obtain after some algebra a single poly- 
nomial equation in single variable x. Thus the prob- 
lem of finding optimal Gaussian measurement on a sin- 
gle mode of a generic three-mode mixed Gaussian state 
localizing maximum entanglement between the remain- 
ing two modes boils down to the finding of roots of a 
high-order polynomial in a single variable which can be 
solved efficiently on a computer. Note, that besides the 
roots (|4"7)) one also has to consider the root y = that 
corresponds to projection of mode C onto vacuum and 
the boundary case y = ±1 corresponding to homodyne 
detection on mode C. These measurements can also po- 
tentially represent the optimal entanglement localization 
strategy. 



B. Optimal homodyne detection 

Homodyne detection can be obtained as a limiting case 
when r — + 00. In this limit Eqs. (j3"6")l and (|4.jp simplify to 



s (ho 
A (ho 



= detA + detB - 2detL> 
det7ABcTc,n(0) 



X22(9) 



(48) 
(49) 



where x (0) = U T (9) X U{9) and r c (0) = U T {9)T C U{9). 
The extremal equation df /d9 = leads to the equation 



(|44[). where j — 9, which can be further arranged into 
the form 



det7 ABC 
c 3 (l + x 2 f 



(50) 



where x = tan 9 and the coefficients hi can be found in 
the Appendix. Inspection of Eq. (j3U|) reveals that the 
problem of finding optimal phase in homodyne detection 
on mode C that maximizes entanglement between modes 
A and B reduces to the problem of finding roots of only 
fourth-order polynomial which can be solved analytically. 
In addition to the roots one has to consider also the lim- 
iting case x — > 00 corresponding to the phase 9 = it/ 2 as 
a possible candidate for the extremum. 



C. Isotropic states 

Previous results reveal that in the case of a general 
Gaussian measurement the maximum localizable entan- 
glement can be found only numerically. However, several 
exceptions exist such as pure three-mode states [l8| or 
bisymmetric states investigated in Sec. IIV1 where fully 
analytical solution can be derived. Besides these states, 
there is yet another class of the so-called isotropic states 
[29l ] for which the optimization task can be resolved ana- 
lytically. The three-mode isotropic states with CM ~fABC 
are Gaussian states that can be transformed by a suit- 
able symplectic transformation S into the normal form 
S"iABcS T = v\, where v > 1 is the three- fold degenerate 
symplectic eigenvalue. Since v = 1 for pure states, the 
CM of an isotropic state can be expressed as 



1ABC = vip, 



(51) 



where 7 P = S^ 1 (S T ) 1 is CM of a pure state. Making 
use of the relation (^7 P ) 2 = —1 [13] that is valid for pure 
states we get 



7 P 1 = ~^7p^- 



(52) 



Hence, we can calculate the inverse matrix (|39[) for 
isotropic state as 



1 



ABC 



1 



;fljABcQ, 



(53) 



where Eqs. |52|) and ([ST]) were used. Assuming the sub- 
matrix C to be again in the standard form C = cl and 
taking into account the relation ft — ©f =1 </ we finally 
arrive at the submatrix Tc for an isotropic state in the 
form: 



(54) 



Consequently, for isotropic states (with C brought into 
the standard form) the parameters g\ and gi vanish and 
the symplectic invariant (|43p is independent of the phase 
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9. This observation greatly simplifies optimization with 
respect to the phase. Since only the invariant S depends 
on 9 the extremal equation df/d9 — is equivalent 
to the equation dS/d9 = which is obviously solved 
by cos(2 6>) = ±x 2 /||xll and sin(20) = ±xi/||xll, where 
Hxll = VXi+xf- Calculating now the second derivative 
d 2 f/d9 2 one finds it is nonnegative for the plus sign and 
therefore the lower symplectic eigenvalue ji is minimized 
by 

S[n{20) ^w cos{20) ^w (55) 

We see that for isotropic states we can find analytically 
optimal phase 9 of the CM even without resorting to 
the limit of infinitely large r corresponding to homodyne 
detection. In this limit we get in particular 



(horn) 
fhnin 



1 



c(hom) 



fr(homh 5 
V min / 



(56) 



where 



1.6 
1.4 
1.2 



3 

<D 1.0 



0.8 
0.6 
0.4 
0.2 - 
0.0 



after localization 
before localization 
homodyne detection 



1.0 



1.5 



2.0 



2.5 

v 



3.0 



3.5 



4.0 



FIG. 1: Dependence of entanglement of modes A and B on u 
for the three-mode Gaussian isotropic state with generating 
covariance matrix (|59[l . The logarithmic negativity of the two- 
mode state before (dotted curve) and after (solid curve) the 
localization is plotted. The shaded area indicates the region 
where homodyne detection (dashed curve) on C is the optimal 
entanglement localization strategy. Otherwise, projection of 
mode C onto coherent states is optimal. 



S^=detA. 



detB 



2detZ? + -(Trx+ llxll), (57) 
2c 



and Eqs. (gU), (HJ) and flM]) were used. 

In the case of projection onto the state with finite 
squeezing r we have to solve the extremal equation 
df jdr = which leads to the equation (03]) with j = r. 
Substituting from Eqs. ([5*5*)) and ([5*1)) into Eqs. (["**"")) and 
(|43p and inserting the obtained symplectic invariants into 
the extremal equation we finally arrive at the polynomial 
equation for optimal squeezing in the form 



(58) 



i=0 



where y — tanhr and the coefficients qi can be found 
in the Appendix. Thus, the problem of finding opti- 
mal Gaussian measurement localizing maximum entan- 
glement for a generic three-mode isotropic state requires 
to calculate roots of only fourth-order polynomial which 
can be performed analytically. 

The isotropic states are of particular interest because 
they represent examples of states for which it is not in 
general optimal to perform homodyne detection. Instead, 
there are isotropic states for which maximum entangle- 
ment is localized by projection onto a finitely squeezed 
state. Consider the isotropic state with CM (["""I)) , where 
v = 2 and 



7p 



/ 3 2 2 \ 
3 0-2 0-2 
2 2 1 
0-2 2 1 
2 1 2 

\ -2 1 2 / 



(59) 



Before measurement on mode C the logarithmic nega- 
tivity of the reduced bipartite state of modes A and B 
reads = 0.189 e-bits. Further, for the CM (J59)) we get 
Hxll = and Eqs. ([5"6"|) and (f5T|) give a lower eigenvalue 
^2 min = 0-636 corresponding to the optimal homodyne 
detection, yielding i?L,hom = 0.654 e-bits. By analyti- 
cally solving polynomial equation (|58p we find that it is 
optimal to take r = which corresponds to projection 
onto the coherent state for which we get /i2, m in — 

0.592 

and El q — 0.757 e-bits. The largest entanglement is 
thus localized by projection onto the coherent state with 
CM 7<¥ = 1 which can be implemented by eight-port 
homodyne detection. For illustration, the dependence of 
the localizable entanglement on the parameter v is plot- 
ted in Fig. CD 



VI. NON-GAUSSIAN MEASUREMENTS 

So far we have limited ourselves to Gaussian measure- 
ments. Recently, it was shown in [l8| that this may 
not be optimal strategy if we allow also non-Gaussian 
measurements. This was illustrated on the example of 
a three-mode pure Gaussian state prepared by mixing a 
vacuum state in mode A with mode B of a two-mode 
squeezed vacuum state 



\X) BC = Vl-A 2 ^A>,n>BC 



(60) 



of modes B and C (see Fig [2]), where A is the squeezing 
parameter. Then it was demonstrated that by perform- 
ing ideal photon number measurement on mode C one 
can localize a higher entanglement between modes A and 
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SPD 




FIG. 2: Scheme demonstrating superiority of non-Gaussian 
localization of entanglement over optimal Gaussian localiza- 
tion. BSas: balanced beam splitter; BScd: unbalanced 
beam splitter with transmittance 77: SPD: ideal single-photon 
detector. See text for details. 

B than wc would get by using optimal Gaussian measure- 
ment, i.e. homodyne detection [lEjj . While the photon 
number detection is suitable for the proof-of-principle of 
superiority of non-Gaussian measurements it is more re- 
alistic to consider the imperfect single-photon detector as 
depicted in Fig [2] Here we show that this non-Gaussian 
measurement strategy still outperforms (in some region 
of the parameter A) the homodyne detection. 

An ideal single-photon detector (SPD) has two out- 
comes, either a click described by the projector 111 = 
1 — |0)(0|, or no click described by the projector onto the 
vacuum state IIo = |0)(0|. Imperfect SPD with detec- 
tor efficiency 77 is then modeled by an unbalanced beam- 
splitter with transmittance rj followed by the ideal de- 
tector (see Fig. [5]). No click of the detector occurs with 
probability po and heralds preparation of modes A and 
B in the mixed state po whereas click is detected with 
probability p\ = 1 — po and prepares the two modes in 
the state p\. 

We quantify the entanglement thus localized between 
modes A and B by the average quantity 

El,ng = PoEsf[po] +PxEm\px}. (61) 

The logarithmic negativity of an arbitrary generally non- 
Gaussian state can be numerically calculated as Ej^[p] = 
log 2 (1 + 2 1 Yli e i|)i where the sum goes over all negative 
eigenvalues of the partially transposed matrix p TA . 
Projection of the mode C onto the vacuum state pre- 
pares the modes B and D in a two-mode squeezed vac- 
uum state that gives after tracing over the auxiliary mode 
D the mode B in a thermal state. Mixing of such a state 
with the vacuum state |0)a on a BSab creates the state 
Po that is obviously separable and hence E_\f[po] = 0. It 
therefore suffices to calculate the logarithmic negativity 
for pi. Describing the beam splitters BSab and BScd 
by the unitary operators Uab and Ucd, respectively, the 
initial state \ipi n ) = |0)a|A)bc|0)d is transformed to the 
state \iP)abcd — UabUcd\^P\i\)- The sought condition- 
ally prepared (unnormalized) state then can be calcu- 



0.4 




FIG. 3: Gaussian localizable entanglement El,g and average 
entanglement El,ng localized by means of the single-photon 
detector with detector efficiency 77 are plotted versus the pa- 
rameter A. See text for details. 

latcd as 

Pi = Tr C B [\iP)abcd(iP\ ( Aab ® H ljC ® Id)] , (62) 

and the corresponding probability reads p\ = Tr ab (pi)- 
After some algebra we get explicitly the normalized den- 
sity matrix p\ — p\/p\ to be 

(l — X 2 ) 00 

Pi = A 1 £ A 2 " [1 - (1 - V 2 r] IVn/AS^U 

(63) 

where 



\i>n) A B = ^J2^ k )\k,n-k) AB , (64) 

\2 2 

and pi = i_^2(i_^2) ■ The negative eigenvalues of the 

matrix p\ A with elements ( p\ A ) = (pi) k a a are cal- 

culated numerically from Eq. (|63[) and the obtained av- 
erage logarithmic negativity (16 lj) is depicted in Fig. [3] 
For comparison we also display in the figure maximum 
Gaussian localizable entanglement [l8[ that is given by 
the entropy of entanglement 

El,g = («a + 1) log 2 («A + 1) - n A log 2 (ru), (65) 
where 

nA = \fr-^- l r (66) 

and that is achieved by homodyne detection. The av- 
erage logarithmic negativity E^^q is plotted in Fig. [3] 
as a function of the squeezing parameter A for various 
values of detector efficiency 77. We can see that even for 
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t] < 1 there is a region of squeezing parameters A where 
El.ng > E LiG so the imperfect SPD still outperforms 
homodyne detection. The region of better performance, 
however, gets smaller with decreasing detector efficiency 
r\ and the maximum value of the difference El,ng~ ^£,G 
reduces and shifts towards smaller values of A. The ad- 
vantage of the SPD over homodyne detection can be best 
understood in the limit of low squeezing A <§; 1 . It follows 
from Eq. (66]) that n A ~ A 4 /4 and E L . G = 0(A 4 ). The 
SPD prepares with probability pi w 2ij 2 X 2 almost pure 
singlet state |V'i)ab = (|01)ab + |10)ab)/v / 2 containing 
one e-bit of entanglement between A and B and therefore 
E l ,ng = 0(\ 2 ). 

VII. CONCLUSIONS 

In the present paper we have analyzed in detail the lo- 
calization of entanglement of mixed multimodc multipar- 
tite Gaussian states by local Gaussian measurements and 
classical communication. In particular, we have found 
that, in contrast to pure Gaussian states, the optimal 
Gaussian measurement that localizes maximum entan- 
glement need not be homodyne detection (i.e. projection 
onto infinitely squeezed states) but rather projection on a 
finitely squeezed state. It is important to stress that any 
such measurement is experimentally feasible with present 
technology and can be implemented on light by splitting 
the signal beam on an unbalanced beam splitter and us- 
ing two balanced homodyne detectors to measure the am- 
plitude quadrature of the first output mode and the phase 
quadrature of the second output mode, respectively. 

The determination of the optimal Gaussian measure- 
ment for localization of entanglement of mixed states is a 
non-trivial problem. Nevertheless, we were able to obtain 
fully analytical results for generic mixed fully symmetric 
Gaussian states, where we proved that homodyne detec- 
tion of the amplitude or phase quadrature is optimal. For 
three-mode mixed Gaussian states we have shown that 
the optimization problem reduces to calculation of roots 
of certain polynomials, which can be very efficiently per- 
formed numerically. The problem simplifies considerably 
if we restrict ourselves to localization by homodyne detec- 
tion where the optimal phase of the homodyne detection 
can be found as a root of only fourth-order polynomial 
which can be calculated analytically. Fully analytical re- 
sults can be obtained also for the three-mode isotropic 
Gaussian states because in this case the determination 
of Gaussian localizable entanglement again boils down 
to finding roots of the fourth-order polynomial. 

Our results contribute to the understanding of the 
structure and properties of entanglement of multimodc 
Gaussian states and we hope that our work will stimulate 



the experimental investigations of entanglement localiza- 
tion for continuous- variable systems. 
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APPENDIX: POLYNOMIAL COEFFICIENTS 

The coefficients of the polynomial ([50)1 read as 

hi = T c ,22Xi + Xi9iu + cg 2 det-f A BC, 

h 3 = 2x\g%u + 2x29\u + 4r c ,22XiX2 + ^cg^det-fABC, 

h 2 = 4r c ,22X2 - Xi Trr c + 4X232U - Xi5i ( u + v ) 

+2c {2gl - gj) det lABC , 
In = -4r c ,nXiX2 - 2i> (X291 + X1S2) - 4:cg 1 g 2 detjABC, 
ho = r c ,nXi + Xi9iv + cgldet^ABc, 

where 

u = xn + el, v = X22 + cl, 

and 

I = detA + dctB - 2detL>. 
The coefficients of the polynomial ([55)) read 
q 4 = \\ X \\ 2 b-a-, 

93 = \\x\\[2a-b-d+ - (a+b- + a-b+)d-], 
q 2 = (a-b+-a+b-) 2 -||x|| 2 (a-&++a+M 

+ (cLa+ — d+a_)(rf_6 + — rf + &_), 
qi = \\x\\[2a + b + d^ - (a + b- + a-b + )d + ], 
9o = ||x|| 2 &+a+, 

where 

a± = ±v 2 {v 2 ±c) 2 1 
b ± = ±(l±c) 2 , 

d± = ±(l±c) 2 /±(det£-detF) 2 -|-(l±c)Trx. 
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